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*vj i In wireless networks, the knowledge of nodal distances is essential for several areas such as system 



configuration, performance analysis and protocol design. In order to evaluate distance distributions in 
random networks, the underlying nodal arrangement is almost universally taken to be an infinite Poisson 
point process. While this assumption is valid in some cases, there are also certain impracticalities to this 
model. For example, practical networks are non-stationary, and the number of nodes in disjoint areas 
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j^ ' are not independent. This paper considers a more realistic network model where a finite number of 

nodes are uniformly randomly distributed in a general d-dimensional ball of radius R and characterizes 
the distribution of Euclidean distances in the system. The key result is that the probability density 
function of the distance from the center of the network to its n* nearest neighbor follows a generalized 
beta distribution. This finding is applied to study network characteristics such as energy consumption, 
interference, outage and connectivity. 
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I. Introduction 
A. Motivation 

In wireless channels, the received signal strength falls off with distance according to a power 
law, at a rate termed the large scale path loss exponent (PLE) [1]. Given a transceiver distance 
d, the signal power at the receiver is attenuated by a factor of d~", where a is the PLE. 
Consequently, in wireless networks, distances (or equivalently, the distribution of distances 
for random networks) between nodes strongly impact the signal-to-noise-and-interference ratios 
(SINRs), and therefore the link reliabilities as well. The knowledge of the nodal distances is 
therefore essential for several important areas such as the throughput and performance analysis 
and the design of protocols and algorithms. 

In wireless networks, nodes are typically scattered randomly over an area or volume, and the 
distance distributions follow from the spatial stochastic process governing the locations of the 
nodes. For the sake of analytical convenience, the arrangement of nodes in a random network 
is commonly taken to be a homogeneous (or stationary) Poisson point process (PPP). For the 
so-called "Poisson network" of intensity A, the number of nodes in any given set of Lebesgue 
measure V is Poisson with mean XV, and the number of nodes in disjoint sets are independent 
of each other. Even though the PPP assumption can lead to some insightful results, practical 
networks differ from Poisson networks in certain aspects. First, networks are formed by usually 
scattering a fixed (and finite) number of nodes in a given area (or very close to it). If nodes are 
uniformly randomly distributed, they form a binomial point process (BPP), which we describe 
shortly. Secondly, since the area or volume of deployment is necessarily finite, the point process 
formed is non-stationary and often non-isotropic, meaning that the network characteristics as 
seen from a node's perspective is not homogeneous for all nodes. For example, the mean farthest 
neighbor distances are larger for nodes lying near the boundary than for the ones closer to the 
center. Furthermore, the number of nodes in disjoint sets are not independent but governed by a 
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multinomial distribution. Fig. \T\ shows a realization of the two processes with the same density 
and depicts that the PPP is clearly not a good model, since there may be more points in the 
realization than the number of nodes deployed. In particular when the number of nodes is small, 
the Poisson model is inaccurate. The main shortcoming of the Poisson assumption is, however, 
the independence of the number of nodes in disjoint areas. For example, if all users or nodes 
are located in a certain part of the network area, the remaining area is necessarily empty. This 
simple fact is not captured by the Poisson model. This motivates the need to study and accurately 
characterize finite uniformly random networks, in an attempt to extend the plethora of results 
for the PPP to the often more realistic case of the BPP We call this new model a "binomial 
network". 




Fig. 1. (Left) A realization of 10 sensor nodes uniformly randomly distributed in a circular area of unit radius. (Right) This 
realization of the Poisson network with the same density (A =3.18) has 14 nodes. The shaded box at the origin represents the 
base station. 



Formally speaking, a rf-dimensional BPP is formed as a result of distributing A^ points 
independently uniformly in a compact set W CM.'^ [3]. The intensity of this process is defined to 
be N/udiW). For any Borel subset V of W, the number of points in V, ^(V), is binomial(n, p) 
with parameters n = N and p = i'd{y)/^d{W), where Vd{) is the standard d-dimensional 
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Lebesgue measure. Accordingly, 

Pr($(y) = k)= ( ', J/(l-p)' 

By this property, the number of nodes in disjoint sets are joint via a multinomial distribution. 
This network model applies to mobile ad hoc networks and wireless networks with infrastructure, 
such as cellular telephony networks or sensor networks. 

In this paper, we study the Euclidean distance properties in a general (i-dimensional bino- 
mial network. We also provide results on network characteristics such as energy consumption, 
interference, outage and connectivity based on the distance distributions. 

B. Related Work 

Even though the knowledge of node locations in wireless networks is important for their 
analysis and design, relatively few results are available in the literature in this area. Moreover, 
much of existing work deals only with moments of the distances (means and variances) or 
characterizes the exact distribution only for very specific system models. 

In [4], the probability density function (pdf) and cumulative distribution function (cdf) of dis- 
tances between nodes are derived for networks with uniformly random and Gaussian distributed 
nodes over a rectangular area. [5] studies mean internodal distance properties for several kinds 
of multihop systems such as ring networks, Manhattan street networks, hypercubes and shuf- 
flenets. [6] provides closed- form expressions for the distributions in (i-dimensional homogeneous 
PPPs and describes several applications of the results for large networks. [7] derives the joint 
distribution of distances of nodes from a common reference point for planar networks with a 
finite number of nodes randomly distributed on a square. [8] considers square random networks 
and determines the pdf and cdf of nearest neighbor and internodal distances. [9] investigates 
one-dimensional multihop networks with randomly located nodes and analyzes the distributions 
of single-hop and multiple-hop distances. 
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II. Distance Distribution to the Neighbors 

In this section, we determine the distribution of the Euclidean distance to the n"^ nearest point 
from the origin for a general rf-dimensional isotropic BPP. It will be established that this random 
variable follows a generalized beta distribution. We also comment on the distances to the nearest 
and farthest neighbors and the void probabilities. 

Theorem 2.1: In a point process consisting of A^ nodes uniformly randomly distributed in a 
(i-dimensional ball of radius R centered at the origin, the Euclidean distance Rn from the origin 
to its ra* nearest neighbor is distributed as a generalized beta distribution i.e., 

, ,, dB{n-l/d+l,N-n + l)^f/r\<i 1 , ,, A r^^. 

where P{x;a,b) is the beta density functioiu defined as P(x;a,b) = ^/ ^s X°~^(l — x)*~^. 
B{a, h) denotes the beta function and is expressible in terms of gamma functions as -B(a, h) = 
T{a)T{b)/T{a + b). 

Proof: Consider the BPP with N points uniformly randomly distributed in a d-dimensional 
ball W of radius R centered at the origin o, i.e., W = Md{o, R). The volume of this ball Vd{W) 
is equal to Qi?'^ [3], where 

Cd 



r(i + t//2) 

is the volume of the unit ball in W'-. The important cases include Ci = 2, C2 = tt and C3 = 47r/3. 
The density of this process is equal to N/{cdR'^). 

The complementary cumulative distribution function (ccdf) of Rn is the probability that there 
are less than n points in the ball 18^(0, r) : 

^«.W = E( JP'(1-P)'^~'' 0<r<i?, (1) 

'Mathematica; PDF[BetaDistribution[a, b],x]. 
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where p = Cdr'^/cdR'^ = (r/R) . 

Fr^ can be written in terms of the regularized incomplete beta function as 

Ffi„(r) = Ji_p(iV-n + l,n), 0<r<i?, (2) 

where 

The pdf of the distance function is Jr^ = — dF^^/dr and we have 



d{l-p)\ (l-p)^->^ 



N~n^n~l 



dr J B{N -n + l,n) 



R\RJ B{N -n + l,n) 

d{l- p)N~npn~l/d 

RB{N -n + l,n) 

dB(n-l/d+l,A'-n+l) /.ry,^_l ^_^_^^ (3) 



R B{N ~n + l,n) \\RJ ' d 

for < r < R. The final equality casts i?„ as a generalized beta-distributed variable. D 

Corollary 2.2: For the practical cases of d = 1 and d = 2, we have 



/i?„(r) = i/3(^;n,iV-n + l) 



i?' \R 
and 

respectively. 

Fig. [2l plots the distance pdfs for the cases of d = 1 and d = 2. 
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Remarks: 

1) The void probability p^ of the point process is defined as the probability of there being 
no point of the process in the test set B C W [3]. For a BPP with N points distributed 
over a set W, it is easy to see that 

p% = il-i.iB)/uiW)f. (4) 

For the isotropic BPP considered above, when the test set is B = Brf(0,r), we have 

p% = h^,iN,i) = ii-pr. 

2) Of interest in particular are the nearest and farthest neighbor distances. The nearest neigh- 
bor distance pdf is given by 



and the distance to the farthest point from the origin is distributed as 

Both these are distributed as the generalized Kumaraswamy distributions [10]. 
3) For a one-dimensional BPP, //j„(r) = /r^_„^i(-R — r), and therefore knowledge of the 
distance pdfs for the nearest \N/2] neighbors gives complete information on the distance 
distributions to the other nodes. 
We wish to compare the distance distributions for a BPP and a PPP. However, note that in 
general, the PPP may have fewer points than the number dropped. In order to make a fair 
comparison, we condition on the fact that there are at least N points present in the PPP model. 
The following corollary establishes the distance pdfs for such a conditioned PPP. Also note that 
conditioned on there being exactly N points present, the PPP is equivalent to a BPP [3]. 
Corollary 2.3: Consider a PPP of density A over a finite volume Md{o,R). Conditioned on 
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there being at least N points in the ball, the distance distribution from the origin to the n* 
nearest neighbor (n < N) is given by 

/rW = ^-w-^ ^-^ —^ re[0,R], (7) 

where Ak{r) := e"^^'*''' {Xcy^ /k\ and fifc(r) := e-^'^'i^'-'"') (Aq {R'^ - r'^))'" /k\. 
Proof: The complimentary conditional cdf of Rn is given by 

F'«„(r) = Pr(<l>(Brf(o,r))<n|<l>(Brf(o,i?))>iV) 
Pr ($ (Brf(o, r)) < n, $ (Brf(o, R)) > N) 



Pr(<l>(Brf(o,i?))>iV) 

(a) 



E::; Pr ($ (B,(o, r)) = A;) Pr ($ (B,(o, i?) \ B,(o, r)) > N - k) 



Pr($(B,(o,i?))>iV) 



where (a) is obtained from the property that the number of points in disjoint sets are independent 
of each other. It is easy to see that 



dr 

and 

dr 

Therefore, we have 

dr 



Ak{r) = { (i) 

-Xdcy~^Ao{r) k = 



d _ Xdcy-HBi{r) - Bi^,{r)) l>0 

—Bi{r) = { (u) 

Xdcdr'^"^Bo{r) I = 0. 



, N-k-l 

- J2 Bi{r) = Xdcy~'BM-k-i{r). (iii) 

=0 



The details of the remainder of the proof are straightforward but tedious and are omitted here. 
Since the pdf of the conditional distance distribution is /jj^ = —dF'j^^Jdr, one basically has to 
differentiate the numerator in (O, and after some simplifications using (i)-(iii), it will be seen 
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that the conditional distance pdf is identical to (|7]). 



D 



Fig. |2] depicts the conditional pdfs of the distances for one- and two-dimensional Poisson 
processes and compares it with the distance distributions for a BPP. 



N = 10, R = 10, d = 1 



N = 10, R = 10,d = 2 



BPP 
Conditioned RPR 





Fig. 2. Distance pdfs for each of thie neighbors for one- and two- dimensional binomial and conditioned Poisson networks. 



When a large number of nodes are distributed randomly over a large area, their arrangement can 
be well approximated by an infinite homogeneous PPP The PPP model for the nodal distribution 
is ubiquitously used for wireless networks and may be justified by claiming that nodes are 
dropped from an aircraft in large numbers; for mobile ad hoc networks, it may be argued that 
terminals move independently of each other. We now present the following corollary to the 
earlier theorem, which also reproduces a result from [6]. 

Corollary 2.4: In an infinite PPP with intensity A on W'-, the distance i?„, between a point 
and its n^^ neighbor is distributed according to the generalized Gamma distribution. 



fR^{r) 



rVin) 



r G 



(9) 



April 29, 2008 



DRAFT 



10 



Proof: If the total number of points N tends to infinity in such a way that the density 
A = N/^CdR^) remains a constant, then the BPP asymptotically (as R -^ oo) behaves as a PPP 



[3]. Taking R = ^N/cdX and applying the limit as A^ ^ oo, we obtain for a PPP, 



fRnir) = lim 



N-.00R T{N -n + l)T{n) 



■{XcdrT lim 1- 



rT(n) ' ' N^oo \ N iV" 



N 



d ,,_d,nr,^(, Xcdr'^y N{N-l)...{N-n + l) 



e 



,,,^,r^ d{\cyr 



rT{n) 

for r G M. This is an alternate proof to the one provided in [6]. D 

For a PPP, we can also specify how the intemodal distances are distributed. By the stationarity 
of the PPP and Slivnyak's theorem [3], the distance between the origin and its n* nearest neighbor 
is the same as the distance between an arbitrary point and its ra* nearest neighbor. [6] studies 
distance distributions in a PPP in more detail and shows that the nearest neighbor distribution for 
any node is exponential and Rayleigh-distributed for the one-dimensional and two-dimensional 
PPP respectively. However, since the BPP is non- stationary, it is not straightforward to specify 
how the internodal distances are distributed in that case. 



April 29, 2008 DRAFT 



11 



III. Moments of R^ 

We now use the distance pdf to compute its moments. The 7* moment of Rn is calculated 
as follows 



nm 



d 



R 



(a) 



R^ r^ 

B{N-n + l,n)Jo 

R^ 
B{N -n + l,n] 



R 



nd—l 



^ \ d 

R 



N-n 



dr. 



r+T/'^"^(l-t)^-"dt 
B,{n + ^/d,N-n + l)\l 



(b) 



R^r{N+l)r{-f/d+n) 
r(n)r{-y/d+N+l) 

00 



if ra + 7/d > 
otherwise 



(10) 



i?7^[7/d]/(iv+l)[T/d] ifn + 7/rf>0 
00 otherwise, 

where Bx[a,b] is the incomplete beta functiorO and xt"' = r(x + n)/T{x) denotes the rising 
Pochhammer symbol notation. Here, (a) is obtained on making the substitution r = Rt^/'^ and 
(h) using the following identities: 



5o(a, h) 



7^e(a) > 

-cx) 7^e(a) < 0, 



and Si(a, h) = B{a, h) if 7^e(6) > 0. 

The expected distance to the n* nearest neighbor is thus 



E(K) 



Rn^^l'^ 
(Ar + i)[i/d] 



(11) 



^Note that 7 £ R in general, and is not restricted to being an integer. 
Mathematica; Beta[x, a, b]. 
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and the variance of R^ is easily calculated as 

Remarks: 

1) For one-dimensional networks, E[_R„] = Rn/{N + 1). Thus, on an average, it is as if the 
points are arranged regular lattice. In particular, when N is odd, the middle point is located 
exactly at the center on average. 

2) On the other hand, as rf — i> oo, E[i?„] -^ R and it is as if all the points are equidistant at 
maximum distance R from the origin. 

3) In the general case, the mean distance of the n^^ neighbor varies as n^^'^' for large n. This 
follows from the series expansion of the Pochhammer sequence [11] 

nM = n'? (1 - C (1/n)) . 

Also, the variance ^ as n increases for d > 2. This is also observed in the case of a 
Poisson network [6]. 

4) The mean internodal distance between the z* and j* nearest neighbors from the origin is 
simply given by (assuming i < j) 

5) For the special case of j/d G Z, we obtain 

IV. Applications to wireless networks 

We now apply the results obtained in the previous section to binomial networks. For the 
system model, we assume a (i-dimensional network over a ball Md{o,R), where A^ nodes are 
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uniformly randomly distributed. Nodes are assumed to communicate with a base station (BS) 
positioned at the origin. The attenuation in the channel is modeled by the large scale path loss 
function g with PLE a, i.e., g{x) = ||a;||~". The channel access scheme is taken to be the slotted 
ALOHA with contention parameter p. 

A. Energy Consumption 

Since battery life in wireless nodes are limited, energy consumption is a key design issue. 
The energy that is required to successfully deliver a packet over a distance r in a medium with 
PLE a is proportional to r". Therefore, the average energy required to deliver a packet from 
the n* nearest neighbor to the BS is given by (flOl) . with 7 = 0;. This approximately scales as 
7i°/<^ when the routing is taken over single hops. When a < d, its more energy-efficient to use 
longer hops than when the PLE is greater than the number of dimensions. 

B. Localization 

In this section, we investigate conditional distance distributions and study their usefulness to 
localization algorithms, where it is beneficial to use as few beacons as possible for the estimation 
process. 

Assume that the nodes in the network can talk to each other and determine how many other 
nodes are closer to the BS than they are, based on the strength of the received signal from the 
BS. However, the channel may be time-varying and hence, it might not be possible for all the 
nodes to determine their exact locations. In such situations, a few beacons can be employed 
that estimate or even precisely measure their distances from the BS. What can be said about the 
distance distributions of the other nodes given this information? 

Suppose we know that the /c* nearest neighbor is at distance s from the center. Then, clearly, 
the first n — 1 nodes are uniformly randomly distributed in Md{o, s) while the more distant nodes 
are uniformly randomly distributed in Md{o, R)\Ed{o, s). Following ([3]), the distance distributions 
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of the first k — 1 nearest neighbors from the origin can be written as 

fR„{r\Rk = s) = ^ — —j ^ -/^{[-j ■,n-- + l,k-n], n<k 

s B[k — n,n) \\s/ a 

for < r < s, which again follows a generalized beta distribution. 
For the remaining nodes i.e., for n > k, we have in r G [s, R], 



fRAr\Rk = s) = -—h^g{N-n + l,n-k) 



d 
dr 



R^ - s'^ B{N - n + l,n - k) 

where q = (r^ — s^)/{R'^ — s^). 

The moments of i?„ are also straightforward to obtain. Following (flOl) . we see that for n < k 

and n + a/d > 0, 

ga [a/d\ 

nK\Rk = s] = j^—^. (13) 

For n > k, we have 

= Fi [ n — k;n — N, — j]n — k + 1;1,1 j 



{n-k)B{N-n + l,n-k) "\ ' ' d 

where Fi[a; 6i, 62; c; x, y] is the Appell hypergeometric function of two variable^. 

Often, it is easiest to measure the nearest neighbor distance. Give this distance as s, we have 
for n > 1, 

^ , , „ , dr'^-^ i 1 - i W^ } } [ W^ 

fRAnRi = s) 






Rd-s'^ B{N -n + l,n-l) 



Mathematica: AppellFl[a, bl, hi, c, x, y]. 
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for r E [s,R]. Also, the mean conditional distances to the remaining neighbors are 

E[R^\R, = s] = ' FJn-l;n-N,-ln-l;l,l-^ 

[n — l)B[i\ — n + l,n — 1) \ a s"- 

C. Interference 

In order to accurately determine network parameters such as outage, throughput or transmission 
capacity, the interference in the system / needs to be known. The mean interference as seen at 
the center of the network is given by 



/i/ = E 



N 



Y.(TnRnn 
n=l 

pJ2^[R--], 



N 



^E[T„]E[/2;°], 



71=1 



n=l 

where T„ G {0, 1} are i.i.d Bernoulli variables (with parameter p) representing whether the n"^ 
nearest node transmits or not in a particular time slot. 

Setting 7 = —a in (flOl) . we can conclude that the mean interference is infinite for rf < a. 
When the number of dimensions is greater than the PLE, we have 

_ pR~''T{N+l) ^ T{n-a/d) 
^'~r{N + l-a/d)^^ Tin) " 

One can inductively verify that 

V^ r(n - a/d) T{k - a/d) k - a/d ^, ^ „ ,,., 

and we obtain after some simplifications. 



-a 



Hi = — ; , d > a. (15) 

d — a 

The unboundedness of the mean interference at practical values of d and a (i.e., d < a) actually 
occurs due to the fact that the path loss model we employ breaks down for very small distances. 
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In fact, it exhibits a singularity at x = 0. In order to overcome this, some authors employ a 
modified path loss law given by min{l, ||a;||~"}. Employing the modified path loss law, the mean 
interference is given by [12] (assuming R > I) 



Hi 



Npd 
Npd 



1 I fl''-"- 
d """ d-a 



d ^ a 



(16) 



^[l/d + \n{R)] d 



a 



D. Connectivity 

We now study the connectivity properties of the binomial network. Define a node to be 
connected to the origin if the SINR at the BS is greater than a threshold 0. Let the nodes 
transmit at unit power and assume noise to be AWGN with variance A'q. In the absence of 
interference, the probability that the BS is connected to its n"^ nearest neighbor is 

= Pr(i?-" > iVoG) 

= 1 - Pr(i?„ > (iVoe)^^/") 

l-/i_p/(iV-n + l,n) e>i?-"/iVo 

1 e < i?~"/iVo, 



(17) 



where p' = ({NqQ) /R) ■ Fig. [3] plots the connectivity probability in a two-dimensional 



binomial network with 25 nodes. 



E. Outage Probability 



An outage O is defined to occur if the SINR at the BS is lower than a certain threshold 9. 
Let the desired transmitter be located at unit distance from the origin, transmit at unit power and 
also not be a part of the original point process. Assuming that the system is interference-limited. 
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N = 25, d = 2, a = 4, R = 20, N|^ = 5e-6 




5 10 

Threshold (dB) 



20 



Fig. 3. The probability of the n* nearest neighbor being connected to the BS for a binomial network with 25 nodes. 



we have 



Pr(C) = Pr[l//<0] 
= Pr [/ > 1/6] 



(18) 



Considering only the interference contribution from the nearest neighbor to the origin, a simple 
lower bound is established on the outage probability as 



Pr(C) > Pr (Tii?r° > 1/e) 



p\l 



P 



Qd/o 



N 



<i?" 

e>i?". 



(19) 
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The empirical values of success probabilities and their upper bounds (1191 ) are plotted for different 
values of A^ in Fig. HI As the plot depicts, the bounds are tight for lower values of N and 0, and 
therefore we conclude that the nearest neighbor contributes to a major portion of the network 
interference. However, note that as a decreases, the bound gets looser since the contributions 
from the farther neighbors are also increased. 



R = 25, d = 2, a = 4, p = 0.5 



0.95- 



0.9- 



0.85 



Q- o.e 



=! 0.75 



0.7- 



0.65 



-10 




■ Empirical value 
Upper bound 



5 10 

Threshold (dB) 



Fig. 4. Comparison of exact success probabilities versus their upper bounds for different values of the system parameters. 



F. Other Applications 

We now list a few other areas where knowledge of the distance distributions is useful. 

• Routing: The question of whether to route over smaller or longer hops is an important, 
yet a nontrivial issue [13], [14]. The knowledge of intemodal distances is necessary for 
evaluating the optimum hop distance and maximizing the progress of a packet towards its 
destination. 
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• Path loss exponent estimation: The issue of PLE estimation is a very important and relevant 
problem [15]. Several PLE estimation algorithms are based on received signal strength 
techniques and these require the knowledge of distances between nodes. 

V. Concluding Remarks 

We point out that the Poisson model for nodal distributions in wireless networks is not accurate 
in many practical situations, and consider the more realistic binomial network model. We derive 
exact analytical expressions for the pdfs of the distances to the n"' nearest neighbor from the 
origin in a general (i-dimensional isotropic random network of radius R. We also analytically 
express the moments of these generalized beta-distributed variables in closed-form. Our findings 
have applications in several problems related to wireless networks such as energy consumption, 
connectivity, localization, interference characterization and outage evaluation. 
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